Quantum magnetization plateaux of an anisotropic ferrimagnetic spin chain 
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The magnetization curve of the {S,s) — (1,1/2) ferrimagnetic alternating spin chain with the 
single-ion anisotropy D is investigated with the numerical exact diagonalization of finite clusters and 
size-scaling analyses. The system has a plateau at 1/3 of the saturation moment, which corresponds 
to the spontaneous magnetization for D = 0. Varying D in the 1/3- magnetized ground state under 
the external field along the axis of D, a quantum phase transition is revealed to occur at the 
critical value D/J = 1.114 ± 0.001 where the plateau vanishes. Except for the critical point, the 
plateau is always opening, but the mechanism is different between D < Dc and D > Dc. The 
change of mechanisms is an evidence to clarify that the plateau originates from the quantization of 
magnetization. 

PACS numbers; 75.10.Jm, 75.40. Mg, 75.50.Gg, 75.40. Cx 



I. INTRODUCTION 

The quantization of the magnetization is one 
of interesting phenomena in low-dimensional mag- 
nets. It is detected as a plateau in the magne- 
tization curve. Such plateaux were actually ob- 
served in high-field measurements of several materi- 
als: the S — I bond-alternating chain [Ni2(Medpt)2(//- 
ox)(^-N3)]ClO4-0.5H2O (Medpt=methyl-bis(3- 
aminopropyl) amine) the organic S = 1 ladder 

3 , 3 ' ,5 , 5 '-t etrakis ( A^-t eri-butylaminoxyl) biphenyl, abbre- 
viated BIP-TENO [p^ , the Shastry-Sutherland sys- 
tem SrCu2 (603)2 [pyiO|] and the zigzag double chain 
NH4CUCI3 0Jl^ etc. In addition some theoretical 
and/or numerical analyses predicted that a magneti- 
zation plateau appears in various other systerns; the 
polymerized chains |l3|-|23), the 5 = 3/2 chain 
the frustrated spin ladder p8|-p2t, several generalized 
spin ladders ||3^-|38t], distorted diamond type spin chain 
l^-ig and some layered sytems [E2p3]. Using the 
Lieb-Schultz-Mattis theorem Q|, a general necessary 
condition for the presence of the plateau was derived p4] 



integer , 



(1) 



where S and m are the sum of spins over all sites and the 
magnetization in the unit period, respectively. 

The ferrimagnetic mixed spin chains have lately at- 
tracted a lot of interest among quantum spin sys- 
tems. Recent synthesizing techniques have produced 
a lot of such materials, for example, the bimetallic 
chain MM'(pbaOII)(H2 0)3-nIl2 ^ and the organic 
one {Mn(hfac)2}3(3R)2 Pq] etc. Thus many experimen- 
tal investigations have been done on such ferrimagnets, 
as well as theoretical ones. E^-pSf In the systems two dif- 
ferent spins S and s {S > s) are arranged alternately in 



a line and coupled by the nearest-neighbor antiferromag- 
netic exchange interaction. They has the spontaneous 
magnetization m = 5* — s in the ground state. Since the 
lowest excitation increasing m has an energy gap, the 
magnetization curve has a plateau at m = S* — s [ p4[ . 
Indeed the plateau satisfies the condition of the quan- 
tization of the magnetization (^. The plateau, how- 
ever, is also realized in the classical limit where S and 
s are infinite with the ratio S/s fixed. Thus it is dif- 
ficult to identify the plateau at m = 5 — s as a result 
from the quantization of the magnetization, unlike any 
higher plateaux atm — S — s-\-l, S— s-\-2,---,S + s— 1 
which should not appear in the classical Heisenberg spin 
systems [ ^9|pC[ |. The previous works |^,^ to inves- 
tigate the anisotropy in the exchange interaction re- 
vealed that the quantum effect stabilizes the plateau at 
m = S — s against the XK-like anisotropy. However, 
it is no more than a quantitative difference between the 
quantum and classical systems. In this paper, to show a 
more definite evidence to clarify that the plateau is a re- 
sult from the quantization, we investigate the single-ion 
anisotropy D effect. It is more realistic than the inter- 
action anisotropy even from an experimental viewpoint. 
For example, the recently synthesized {S,s) = (1,1/2) 
chain NiCu(pba)(D20)3-2D20 |p3] would possibly have 
the anisotropy D on the Ni ion JS = 1). Thus we study 
on the (1,1/2) system with D ior S ~ I, using the numer- 
ical diagonalization and the level spectroscopy method 
p5|-|68| under the twisted boundary condition |^,^. It 
will reveal that the system has a quantum phase tran- 
sition with respect to I? at m = 1/2 due to the change 
of plateau formation mechanisms, which clarifies a full 
quantum nature of the plateau. 
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II. CLASSICAL SYSTEM 




The ferrimagnetic {S,s) — (1,1/2) mixed spin chain 
with the single-ion anisotropy under an external mag- 
netic field H is described by the Hamiltonian 



0.5 



H — Ho + TYcx 

N 

n, = Y.{s,-s, + s,-s,+, + D{s]y} 



N 



(2) 
(3) 

(4) 



Throughout this paper, we consider only the case when 
the external field is along the symmetry axis of D. We 
note that D{sj)^ gives a constant D/A. In order to clar- 
ify the quantum nature of the plateau later, we show 
the properties of the classical spin system with the same 
Hamiltonian where S and s are the classical vectors with 
the amplitudes 1 and 1/2, respectively. The magnetiza- 
tion process at T = can be solved with the standard 
variation of the Hamiltonian with respect to the two vec- 
tors S and s. In the isotropic case of the classical sys- 
tem the ground state magnetization curve begins at the 
spontaneous magnetization m ~ S — s — 1/2 and has 
a plateau there, shown as a dot-dashed line in Fig. |^. 
The complete Neel order along H {Sj — l,s| — —1/2) 
is realized in the plateau state, while the canted Neel or- 
der (in the xy plane) occurs for 1/2 < m < 3/2. Thus 
the plateau should originate from the classical Neel or- 
der. Since the Neel order is oriented in the xy plane for 
D > 0, the magnetization curve starts from m = 0. The 
plateau at m — 1/2 still appears for small positive D, 
as a long dashed line {D — 0.05) in Fig. 0, because the 
Neel order along H can be realized there. With increas- 
ing D, however, the plateau disappears at the critical 
value Dc — 0.0572 and it does not appear any more for 
D > Dc, shown in Fig. The breakdown of the plateau 
due to the easy plane anisotropy Z? > is qualitatively 
the same as the case of the XY-\ike coupling anisotropy 
pl| , |62[ . In the quantum system discussed in the following 
sections, however, the plateau will be revealed to appear 
again for larger D due to the quantum effect, in contrast 
to the coupling anisotropy. 
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FIG. 1. Magnetization curves of the classical 

{S,s) — (1,1/2) chain for various D. The plateau appears 
at m = 1/2 only for D < Dc = 0.0572. 



III. QUANTUM MECHANISMS OF PLATEAU 

In the quantum spin system described by the Hamilto- 
nian (|^) it would be efficient to introduce the composite 
spin picture, where S* = 1 is considered as the triplet state 
of two 1/2 spins |7l|. Using this picture we clarify the 
two different mechanisms of the plateau at m — 1/2; the 
Haldane and large- Z? mechanisms for small and large D, 
respectively. The names of these two mechanisms origi- 
nate from the Haldane and large-D phases of the S = 1 
antiferromagnetic chain where the critical point was re- 
vealed to be D ^ 1 |7^-|7^. Similar mechanisms were 
also proposed for the 1/3 plateau of the S = 3/2 chain 

Mm- 



A. Haldane Plateau 

For smaller D each S = 1 site can be in any state of 
=-l, and -1-1. Thus the two 1/2 spins represent- 
ing S = 1 have only to be symmetrized. In the plateau 
state at m = 1/2 (1/3 of the saturation magnetization), 
the antiferromagnetic trimer state, which is schematically 
shown in Fig. g, is expected to be realized. In Fig. || a 
solid ellipse represents the trimer 

i^> = ^(im)-2inT) + im)), (5) 

and a dotted one corresponds to an 5 = 1 site where 
the two 1/2 spins should be symmetrized. When the 
plateau is based on the effective mechanism in Fig. 
we call it the Haldane plateau. In the ideal state where 
this picture is exactly realized, the expectation values of 
the z-component are ^ , — ^ and ^ for the three spins in 
each trimer, respectively. Thus the original system is 
expected to have (5^) = | and (s^) = 



i in the ideal 



2 



Haldane state. It suggests that the Haldane state has the 
Neel order along H , but the amphtude is smaller than the 
classical system ((5^) = 1 and (s^) = — |) because of the 
quantum spin reduction. We note that the D = —oo 
case is not the ideal Haldane case. In this case, classical 
values {S^) — 1 and (s^) = —1/2 will be realized. 




FIG. 2. Schematic picture of the Haldane mechanism of 
the plateau at m = 1/2. Each circle represents a 1/2 spin. 
A solid ellipse is a trimer and a dotted ellipse means sym- 
metrized two spins. 



B. Large-D Plateau 

Each S = 1 site tends to have ^ for larger positive 
D. In the composite spin picture — corresponds to 
one of the triplet states 
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(I n> + UT)) 



(6) 



Thus the large- 1? mechanism of the 1/3 plateau (m = 
1/2) is presented schematically in Fig. ||, where a rect- 
angle represents the triplet state (||) at S* = 1 site and 
each s = 1/2 site has = 1/2. In the ideal large- 1? 
state the expectation values should be {S^) = and 
(s^) = +1/2. Obviously no Neel order is realized along 
z-axis in the large-D phase. The large-D plateau results 
from the quantization of S* = 1. Thus it is never realized 
in the classical spin system. 
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FIG. 3. Schematic picture of the large-D mechanism of 
the plateau at m = 1/2. Each circle represents a 1/2 spin. A 
solid rectangle is one of the triplet (6). 



E(M) is the lowest energy level in the subspace where the 
eigenvalue of ^^■(5'|-fs|) is M and m = 1/2 corresponds 
to M = N/2 for the A^-unit system. A is also the length 
of the plateau. Based on the analogy with the S = 1 
chain, it is expected that the critical point Dc between 
the Haldane and large-D phases is a Gaussian fixed point. 
Thus the gap would vanish just at Dc and open in both 
phases. The scaled gap iVA calculated for several finite 
systems under the periodic boundary condition using the 
numerical diagonalization is plotted versus D in Fig. ^ 
It indicates that there exists a critical point at Z) ^ 1.1 
where the scaled gap A is independent of N, namely 
the system is gapless (A ~ 1/-/V). It also suggests that 
the gap is opening for both sides of Dc, as expected. In 
general, the fixed point of the phenomenological renor- 
malization equation |64[| 



NAn{D) = iN + 2)AN+2iD') 



(8) 



depends on the system size TV for small clusters. How- 
ever, the behavior of the scaled gap in Fig. |^ suggests 
that the fixed point is almost independent of TV, that 
is Dc^ifi ~ Defies- It implies that the higher-order size 
correction {o{l/N)) of the gap A is negligible even for 
N —A, 6, and 8 in the present system. Thus some sys- 
tematic size-scaling analyses would lead to a precise esti- 
mation of the critical point in the thermodynamic limit, 
even based only on such small cluster calculations. 
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FIG. 4. Scaled gap A^'A plotted versus D. It suggests that 
the critical point is around D = 1.1. 



IV. QUANTUM CRITICAL POINT 

The existence of the two different mechanisms of the 
1/3 plateau suggests that there is a quantum phase tran- 
sition between them with respect to the parameter D in 
the ground state at m — 1/2. A useful order parameter 
to investigate the phase transition is the spin excitation 
gap at m = 1/2 

A = E{M + 1)+E{M-1)-2E{M). (7) 



One of the most precise methods to determine the 
Gaussian fixed point of the one-dimensional quantum 
systems is the level spectroscopy ||65|-|6^ under the 
twisted boundary condition ||6^,|7^. The twisted bound- 
ary condition means that the sign of the coupling con- 
stant for the XY component is switched in the exchange 
interaction at the boundary. According to the method, 
the critical point is determined as a crossing point of the 
lowest two energy levels under twisted boundary condi- 
tions. The two levels Ei and E2 are plotted versus D for 
the finite system with =4, 6 and 8 in Fig. ^. The size 
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dependence of the crossing point is so small that a precise 
Dc is expected to be obtained. Since the effect of some 
irrelevant fields [p6[ should yield the size correction pro- 
portional to 1/N^we extrapolate the size-dependent Dc 
to the thermodynamic limit using the plot versus 
in Fig. |. The rcsuh is = 1.114 ± 0.001. The good 
agreement of the size correction with 1/A^^ in Fig. [^jus- 
tifies the accuracy of the present estimation. In fact, the 
plateau-nonplateau phase boundary of some spin ladder 
systems determined by the level spectroscopy even with 
small cluster calculations well agreed with the exact re- 
sult in some ideal limits. ^,|l|,^ The above scaled gap 
analysis also supports the existence of the quantum crit- 
ical point in the present system. 
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FIG. 5. Lowest two energy levels under the twisted bound- 
ary condition. The crossing point is the size-dependent criti- 
cal point. 



for finite systems are plotted versus D in Fig. ^ It in- 
dicates that the sign of (s^) switched around the critical 
point. It means that the Neel order along z-axis exists 
only in the Haldane phase. Thus it is also consistent with 
the schematic pictures of the two plateau mechanisms. 




FIG. 7. Sublattice magnetizations {S^) and (s^) plotted 
versus D. It suggests that the sign of (s^) is changed at the 
critical point. 



V. MAGNETIZATION CURVES 

Finally, we present the ground state magnetization 
curve of the quantum systems described by the Hamilto- 
nian (||) for several values of D, using some size scaling 
techniques [75| ] applied for the numerical energy levels of 
finite systems up to iV = 10. The conformal field theory 
in one-dimensional quantum systems [|7^-|78[| applied for 
the present model predicted that the in gapless phases 
the size dependence of the energy gap have the asymp- 
totic forms 
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FIG. 6. Estimation of the critical point Dc in the thermo- 
dynamic limit based on fitting a line ~ l/N"^ . 

In order to confirm the realization of the two mech- 
anisms of the plateau discussed in the previous section, 
the expectation values of the z component (S*^) and (s^) 



E{M + 1) - E{M) - H{m) + T^v.r,—, 



(9) 



E{M) - E{M -I) H{m) -7:vs7]—, (10) 

where TV and M vary with m = M/N fixed. Thus the 
forms are useful to estimate the magnetic field H for sev- 
eral values of m which can be obtained from all the com- 
binations of N and M for available finite systems. The 
method works except for the plateau. On the other hand, 
at the massive point like plateaux the Shanks transfor- 
mation | |79| , which is a technique to accelerate the con- 
vergence of sequences, is useful. The general form of the 
transformation for a sequence a„ is 



ftn-lfln+l 



-1 + ftn+l — 2a„ 



(11) 



We show the result of the magnetization curve for several 
values of D in Fig. ||, where only the polynomial curves 
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suitably fitted to the obtained points based on the above 
method. Fig. ^ shows that with increasing D up to the 
critical point D = Dc = 1-114: the plateau at m = 1/2 is 
shrinking, while it is opening again for D > Dc, in con- 
trast to the classical system. The plateau is due to the 
Haldane mechanism for D < Dc while the large-D mech- 
anism for D > Dc- This change of mechanisms should 
be a clear evidence of the quantum nature in the plateau 
formation. 



1.5 



0.5 



D=0.0 






D=0.5 


/ / / / 


D=1.114 

D=1.5 


/ / / ' 

/ / / / 

/ / / - 


D=2.0 


/'/ //' / 












) 1 2 


3 4 



H 

FIG. 8. Magnetization curves of the quantum system for 
several values of D. The plateau at m = 1/2 is caused by the 
Haldane mechanism for D = and 0.5, while the large-D one 
for D = 1.5 and 2.0. The curve has no plateau just at the 
critical value D — 1.114. 



VI. SUMMARY 

The magnetization process of the {S, s) = (1, 1/2) ferri- 
magnetic mixed spin chain with the single-ion anisotropy 
D was investigated using the numerical exact diagonal- 
ization and some size-sailing analyses. It revealed that 
the mechanism of the plateau at m — 1/2 is changes 
from the Haldane to large- 1? ones at the Gaussian quan- 
tum critical point Dc = 1.114 ± 0.001. It justifies that 
the plateau at m = 1/2 originates from the quantization 
of the magnetization, although a similar plateau also ap- 
pears in the classical system. 
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